This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.4454148

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.1 while 0.9 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

13

1

1

0.185

0.151

0.035

0.502

1.291

0.126

0.270

NSEI

p_pelagic

13

1

2

0.070

0.051

0.011

0.236

1.370

0.060

0.317

NSEI

p_pelagic

13

10

1

0.229

0.204

0.067

0.523

1.193

0.119

0.210

NSEI

p_pelagic

13

10

2

0.089

0.073

0.020

0.252

1.206

0.062

0.260

NSEI

p_pelagic

13

11

1

0.233

0.209

0.070

0.531

1.177

0.119

0.202

NSEI

p_pelagic

13

11

2

0.091

0.075

0.022

0.247

1.197

0.059

0.271

NSEI

p_yellow

13

11

1

0.407

0.361

0.132

0.965

1.012

0.207

0.222

NSEI

p_yellow

13

11

2

0.409

0.366

0.129

0.969

1.012

0.207

0.208

NSEI

p_pelagic

13

12

1

0.242

0.218

0.072

0.535

1.169

0.121

0.199

NSEI

p_pelagic

13

12

2

0.097

0.081

0.023

0.256

1.196

0.062

0.259

NSEI

p_yellow

13

12

2

0.408

0.362

0.132

0.966

1.011

0.208

0.221

NSEI

p_pelagic

13

13

1

0.254

0.233

0.081

0.558

1.128

0.121

0.173

NSEI

p_pelagic

13

13

2

0.101

0.084

0.027

0.271

1.165

0.064

0.265

NSEI

p_pelagic

13

14

1

0.248

0.230

0.087

0.507

1.112

0.110

0.163

NSEI

p_pelagic

13

14

2

0.102

0.089

0.028

0.249

1.125

0.059

0.220

NSEI

p_pelagic

13

15

1

0.230

0.213

0.079

0.484

1.095

0.105

0.161

NSEI

p_pelagic

13

15

2

0.100

0.087

0.029

0.241

1.116

0.057

0.230

NSEI

p_yellow

13

15

1

0.400

0.353

0.123

0.967

1.012

0.208

0.226

NSEI

p_yellow

13

15

2

0.404

0.358

0.129

0.970

1.011

0.207

0.222

NSEI

p_pelagic

13

16

1

0.269

0.251

0.101

0.523

1.105

0.112

0.161

NSEI

p_pelagic

13

16

2

0.112

0.100

0.033

0.263

1.088

0.059

0.203

NSEI

p_yellow

13

16

2

0.412

0.362

0.132

0.969

1.012

0.208

0.240

NSEI

p_pelagic

13

17

1

0.289

0.273

0.114

0.565

1.077

0.114

0.141

NSEI

p_pelagic

13

17

2

0.121

0.108

0.036

0.284

1.078

0.065

0.199

NSEI

p_pelagic

13

18

1

0.301

0.286

0.121

0.559

1.048

0.114

0.132

NSEI

p_pelagic

13

18

2

0.128

0.114

0.043

0.286

1.070

0.064

0.220

NSEI

p_pelagic

13

19

1

0.311

0.298

0.127

0.570

1.058

0.114

0.114

NSEI

p_pelagic

13

19

2

0.133

0.120

0.044

0.301

1.060

0.066

0.197

NSEI

p_yellow

13

19

2

0.409

0.361

0.121

0.966

1.012

0.208

0.230

NSEI

p_pelagic

13

2

1

0.187

0.154

0.038

0.503

1.326

0.126

0.263

NSEI

p_pelagic

13

2

2

0.073

0.054

0.011

0.248

1.338

0.062

0.306

NSEI

p_pelagic

13

20

1

0.323

0.309

0.138

0.577

1.041

0.113

0.123

NSEI

p_pelagic

13

20

2

0.143

0.129

0.050

0.310

1.044

0.068

0.205

NSEI

p_yellow

13

20

1

0.403

0.355

0.122

0.964

1.013

0.208

0.231

NSEI

p_yellow

13

20

2

0.412

0.367

0.123

0.970

1.012

0.207

0.217

NSEI

p_pelagic

13

21

1

0.347

0.336

0.153

0.600

1.016

0.116

0.095

NSEI

p_pelagic

13

21

2

0.155

0.141

0.052

0.335

1.034

0.073

0.193

NSEI

p_yellow

13

21

1

0.406

0.358

0.126

0.966

1.011

0.206

0.233

NSEI

p_yellow

13

21

2

0.409

0.363

0.132

0.968

1.013

0.208

0.221

NSEI

p_pelagic

13

22

2

0.159

0.149

0.058

0.314

1.024

0.069

0.144

NSEI

p_yellow

13

22

1

0.401

0.356

0.131

0.968

1.011

0.207

0.217

NSEI

p_pelagic

13

23

2

0.107

0.097

0.040

0.221

1.010

0.049

0.206

NSEI

p_yellow

13

23

2

0.459

0.424

0.147

0.975

1.012

0.211

0.166

NSEI

p_yellow

13

26

1

0.364

0.312

0.098

0.958

1.014

0.210

0.247

NSEI

p_yellow

13

27

1

0.356

0.300

0.102

0.956

1.015

0.208

0.269

NSEI

p_yellow

13

28

1

0.360

0.313

0.106

0.940

1.018

0.200

0.235

NSEI

p_yellow

13

28

2

0.451

0.417

0.159

0.967

1.010

0.201

0.169

NSEI

p_pelagic

13

3

1

0.190

0.157

0.037

0.499

1.287

0.124

0.267

NSEI

p_pelagic

13

3

2

0.073

0.056

0.012

0.237

1.284

0.060

0.283

NSEI

p_yellow

13

3

2

0.409

0.366

0.124

0.972

1.011

0.209

0.206

NSEI

p_pelagic

13

4

1

0.196

0.162

0.043

0.521

1.282

0.128

0.266

NSEI

p_pelagic

13

4

2

0.075

0.056

0.014

0.242

1.282

0.061

0.311

NSEI

p_yellow

13

4

2

0.406

0.359

0.127

0.965

1.011

0.209

0.225

NSEI

p_pelagic

13

5

1

0.200

0.168

0.043

0.510

1.260

0.124

0.258

NSEI

p_pelagic

13

5

2

0.076

0.058

0.014

0.249

1.298

0.062

0.292

NSEI

p_yellow

13

5

1

0.403

0.356

0.130

0.963

1.012

0.206

0.228

NSEI

p_pelagic

13

6

1

0.206

0.178

0.047

0.504

1.244

0.121

0.231

NSEI

p_pelagic

13

6

2

0.078

0.060

0.016

0.256

1.266

0.061

0.296

NSEI

p_pelagic

13

7

1

0.211

0.184

0.054

0.509

1.246

0.121

0.223

NSEI

p_pelagic

13

7

2

0.080

0.063

0.017

0.238

1.262

0.059

0.289

NSEI

p_pelagic

13

8

1

0.215

0.188

0.054

0.515

1.270

0.122

0.222

NSEI

p_pelagic

13

8

2

0.084

0.066

0.018

0.256

1.235

0.063

0.288

NSEI

p_pelagic

13

9

1

0.221

0.196

0.058

0.519

1.192

0.119

0.211

NSEI

p_pelagic

13

9

2

0.087

0.069

0.019

0.256

1.235

0.062

0.292

NSEI

p_slope

14

4

1

0.078

0.058

0.011

0.262

1.012

0.068

0.292

NSEO

p_pelagic

16

1

1

0.314

0.294

0.088

0.626

1.146

0.148

0.136

SSEO

p_pelagic

16

1

2

0.296

0.280

0.078

0.602

1.126

0.143

0.112

SSEO

p_pelagic

16

10

1

0.327

0.308

0.106

0.625

1.087

0.141

0.135

SSEO

p_pelagic

16

10

2

0.311

0.291

0.095

0.621

1.100

0.140

0.143

SSEO

p_pelagic

16

11

1

0.330

0.313

0.105

0.641

1.107

0.141

0.120

SSEO

p_pelagic

16

11

2

0.314

0.298

0.098

0.616

1.082

0.137

0.117

SSEO

p_pelagic

16

12

1

0.339

0.325

0.111

0.643

1.081

0.141

0.099

SSEO

p_pelagic

16

12

2

0.317

0.300

0.103

0.633

1.092

0.141

0.121

SSEO

p_pelagic

16

13

1

0.340

0.322

0.118

0.644

1.076

0.140

0.129

SSEO

p_pelagic

16

13

2

0.319

0.307

0.103

0.603

1.091

0.131

0.092

SSEO

p_pelagic

16

14

1

0.326

0.312

0.118

0.617

1.067

0.132

0.106

SSEO

p_pelagic

16

14

2

0.321

0.307

0.109

0.602

1.067

0.132

0.106

SSEO

p_pelagic

16

15

1

0.335

0.324

0.124

0.610

1.064

0.129

0.085

SSEO

p_pelagic

16

15

2

0.325

0.313

0.113

0.617

1.058

0.131

0.092

SSEO

p_pelagic

16

16

1

0.333

0.321

0.128

0.617

1.056

0.128

0.094

SSEO

p_pelagic

16

16

2

0.328

0.315

0.123

0.605

1.042

0.127

0.103

SSEO

p_pelagic

16

17

1

0.391

0.384

0.165

0.657

1.033

0.129

0.054

SSEO

p_pelagic

16

17

2

0.351

0.340

0.128

0.635

1.061

0.134

0.082

SSEO

p_pelagic

16

18

1

0.381

0.369

0.152

0.651

1.037

0.130

0.093

SSEO

p_pelagic

16

18

2

0.351

0.340

0.135

0.622

1.036

0.129

0.085

SSEO

p_pelagic

16

19

1

0.371

0.361

0.150

0.628

1.031

0.127

0.079

SSEO

p_pelagic

16

19

2

0.352

0.340

0.142

0.621

1.031

0.126

0.095

SSEO

p_pelagic

16

2

1

0.317

0.297

0.088

0.641

1.130

0.149

0.134

SSEO

p_pelagic

16

2

2

0.294

0.275

0.078

0.615

1.130

0.144

0.132

SSEO

p_pelagic

16

20

1

0.382

0.376

0.161

0.646

1.012

0.125

0.048

SSEO

p_pelagic

16

20

2

0.364

0.353

0.156

0.632

1.020

0.124

0.089

SSEO

p_pelagic

16

21

1

0.365

0.353

0.163

0.621

1.022

0.120

0.100

SSEO

p_pelagic

16

21

2

0.366

0.355

0.155

0.627

1.016

0.123

0.089

SSEO

p_pelagic

16

22

2

0.379

0.371

0.171

0.633

1.013

0.122

0.066

SSEO

p_pelagic

16

3

1

0.316

0.298

0.091

0.627

1.113

0.147

0.122

SSEO

p_pelagic

16

3

2

0.296

0.275

0.084

0.615

1.127

0.142

0.148

SSEO

p_pelagic

16

4

1

0.320

0.301

0.093

0.632

1.140

0.148

0.129

SSEO

p_pelagic

16

4

2

0.298

0.281

0.084

0.605

1.103

0.140

0.122

SSEO

p_pelagic

16

5

1

0.314

0.298

0.093

0.617

1.136

0.143

0.112

SSEO

p_pelagic

16

5

2

0.299

0.285

0.086

0.603

1.124

0.140

0.100

SSEO

p_pelagic

16

6

1

0.315

0.297

0.091

0.622

1.122

0.146

0.123

SSEO

p_pelagic

16

6

2

0.301

0.279

0.085

0.613

1.114

0.142

0.155

SSEO

p_pelagic

16

7

1

0.324

0.310

0.093

0.641

1.118

0.147

0.095

SSEO

p_pelagic

16

7

2

0.302

0.282

0.092

0.611

1.124

0.140

0.143

SSEO

p_pelagic

16

8

1

0.324

0.307

0.101

0.628

1.122

0.145

0.117

SSEO

p_pelagic

16

8

2

0.303

0.287

0.092

0.607

1.105

0.138

0.116

SSEO

p_pelagic

16

9

1

0.327

0.313

0.101

0.630

1.094

0.143

0.098

SSEO

p_pelagic

16

9

2

0.310

0.292

0.090

0.615

1.093

0.140

0.128

SSEO

p_pelagic

3

1

2

0.191

0.175

0.065

0.391

1.012

0.085

0.189

PWSI

p_pelagic

3

10

2

0.191

0.179

0.066

0.390

1.011

0.084

0.142

PWSI

p_pelagic

3

12

1

0.556

0.563

0.295

0.791

1.017

0.129

-0.054

PWSI

p_pelagic

3

13

1

0.557

0.563

0.296

0.795

1.010

0.131

-0.046

PWSI

p_pelagic

3

13

2

0.191

0.177

0.068

0.387

1.011

0.084

0.167

PWSI

p_pelagic

3

14

1

0.559

0.564

0.305

0.786

1.010

0.127

-0.039

PWSI

p_pelagic

3

15

2

0.188

0.176

0.067

0.392

1.010

0.083

0.144

PWSI

p_pelagic

3

2

1

0.556

0.561

0.297

0.793

1.014

0.128

-0.039

PWSI

p_pelagic

3

4

1

0.555

0.562

0.298

0.787

1.015

0.130

-0.054

PWSI

p_pelagic

3

5

1

0.556

0.559

0.297

0.789

1.011

0.128

-0.023

PWSI

p_pelagic

3

6

1

0.554

0.558

0.295

0.792

1.012

0.132

-0.030

PWSI

p_pelagic

3

6

2

0.192

0.177

0.067

0.396

1.013

0.085

0.175

PWSI

p_pelagic

3

7

1

0.555

0.555

0.293

0.799

1.010

0.132

0.000

PWSI

p_pelagic

3

8

2

0.192

0.179

0.068

0.391

1.011

0.085

0.153

PWSI

p_pelagic

3

9

1

0.555

0.558

0.294

0.790

1.013

0.129

-0.023

PWSI

p_pelagic

4

14

2

0.552

0.559

0.288

0.789

1.010

0.131

-0.054

PWSO

p_pelagic

4

15

2

0.551

0.554

0.292

0.791

1.011

0.129

-0.023

PWSO

p_black

8

40

1

0.491

0.418

0.342

0.964

1.011

0.194

0.376

afognak

p_black

8

48

1

0.536

0.491

0.411

0.971

1.010

0.148

0.303

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta2_pH

11

1

0.425

0.364

0.181

1.053

1.101

0.237

0.257

CSEO

beta2_pelagic

11

2.538

1.641

0.252

7.081

1.505

2.180

0.411

CSEO

beta0_yellow

11

-0.225

-0.186

-0.589

0.002

1.215

0.229

-0.170

CSEO

beta0_pelagic

13

-2.266

-2.141

-4.534

-0.261

1.429

1.194

-0.105

NSEI

beta1_pelagic

13

5.009

4.654

1.598

8.793

1.378

2.167

0.164

NSEI

beta2_pelagic

13

0.119

0.092

0.051

0.316

1.300

0.080

0.340

NSEI

beta3_pelagic

13

29.692

28.694

23.254

41.049

1.206

5.007

0.199

NSEI

beta2_pelagic

14

1.317

0.699

0.221

5.320

1.311

1.437

0.430

NSEO

beta2_yellow

14

-5.908

-4.432

-18.626

-0.379

1.193

4.975

-0.297

NSEO

beta0_pH

15

4

0.838

1.011

-0.477

1.524

1.182

0.563

-0.307

SSEI

beta1_pH

15

4

0.911

0.735

0.220

2.314

1.182

0.603

0.292

SSEI

beta0_pelagic

16

-0.940

-0.864

-2.242

-0.016

1.324

0.634

-0.120

SSEO

beta1_pelagic

16

2.368

2.297

0.857

4.789

1.446

1.037

0.068

SSEO

beta2_pelagic

16

0.421

0.198

0.079

2.725

1.134

0.698

0.319

SSEO

tau_beta2_pH

2

3

0.250

0.111

0.042

1.246

1.168

0.710

0.196

NG

tau_beta1_pH

3

4

2.043

0.085

0.011

16.263

1.103

6.008

0.326

PWSI

beta0_pH

3

1

1.328

1.407

0.395

1.780

1.222

0.337

-0.234

PWSI

beta1_pH

3

1

2.002

1.799

1.130

3.748

1.254

0.669

0.303

PWSI

beta0_pH

8

1

-0.749

-0.663

-2.276

-0.044

1.116

0.492

-0.175

afognak

beta1_pH

9

1

2.134

1.946

1.219

4.436

1.131

0.751

0.250

eastside

lambda_H

10

0.242

0.162

0.027

0.863

1.175

0.493

0.162

northeast

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

614.00

215.00

1,027.00

337.00

Prop. of total

0.74

0.26

0.75

0.25

Prop. median less than 1

0.24

0.26

0.31

0.45

Prop. median less than 10

0.34

0.58

0.50

0.77

Prop. median less than 50

0.48

0.84

0.64

0.84

Median of the median values

60.09

5.49

9.90

1.54

Mean of the median values

581.01

153.15

253.01

168.58

Mean skew of posteriors

0.25

0.29

0.27

0.25

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

54.00

8.00

88.00

29.00

Prop. of total

0.87

0.13

0.75

0.25

Prop. median less than 1

0.04

0.62

0.58

0.69

Prop. median less than 10

0.06

1.00

0.60

1.00

Prop. median less than 50

0.06

1.00

0.68

1.00

Median of the median values

225.37

0.21

0.12

0.18

Mean of the median values

458.82

1.55

120.08

0.87

Mean skew of posteriors

0.27

0.21

0.19

0.18

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_pelagic

2

beta0_pH

3

beta0_yellow

1

beta1_pelagic

2

beta1_pH

3

beta2_pelagic

4

beta2_yellow

1

beta3_pelagic

1

Bp_ay

12

Bp_ayg

9

Bp_ayu

13

eps_pH

1

Hb_ay

13

Hb_ayg

11

Hb_ayu

13

Ho_ay

3

Ho_ayu

5

Hp_ay

15

Hp_ayg

12

Hp_ayu

15

Hs_ayu

1

Hy_ay

3

Hy_ayu

4

lambda_H

1

p_pelagic

43

p_slope

1

R_ay

3

R_ayg

5

R_ayu

5

Rb_ay

4

Rb_ayg

5

Rb_ayu

8

Rd_ay

5

Rd_ayu

7

re_pelagic

5

re_pH

2

Ro_ay

10

Ro_ayg

5

Ro_ayu

30

Rp_ay

3

Rp_ayg

4

Rp_ayu

10

Rs_ayg

1

Rs_ayu

3

Ry_ay

7

Ry_ayg

5

Ry_ayu

24

tau_beta2_pH

1

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.014 0.068 -0.139 -0.018 0.132
mu_bc_H[2] 0.076 0.054 -0.022 0.073 0.188
mu_bc_H[3] -0.271 0.067 -0.399 -0.271 -0.132
mu_bc_H[4] -0.931 0.181 -1.298 -0.925 -0.592
mu_bc_H[5] 1.389 1.272 0.162 1.095 4.607
mu_bc_H[6] -2.403 0.321 -3.013 -2.413 -1.769
mu_bc_H[7] -0.454 0.102 -0.668 -0.451 -0.262
mu_bc_H[8] 0.249 0.352 -0.345 0.217 0.992
mu_bc_H[9] -0.250 0.126 -0.492 -0.250 0.005
mu_bc_H[10] 0.057 0.066 -0.065 0.057 0.192
mu_bc_H[11] -0.124 0.037 -0.196 -0.123 -0.051
mu_bc_H[12] -0.252 0.108 -0.486 -0.247 -0.056
mu_bc_H[13] -0.159 0.073 -0.299 -0.160 -0.014
mu_bc_H[14] -0.309 0.090 -0.493 -0.306 -0.141
mu_bc_H[15] -0.327 0.047 -0.413 -0.328 -0.234
mu_bc_H[16] -0.190 0.366 -0.831 -0.219 0.635
mu_bc_R[1] 1.041 0.120 0.806 1.040 1.273
mu_bc_R[2] 1.016 0.077 0.855 1.020 1.158
mu_bc_R[3] 0.220 0.089 0.041 0.221 0.390
mu_bc_R[4] 0.655 0.198 0.241 0.661 1.019
mu_bc_R[5] 1.682 0.365 0.982 1.674 2.414
mu_bc_R[6] -1.430 0.336 -2.090 -1.431 -0.743
mu_bc_R[7] 0.517 0.166 0.187 0.521 0.832
mu_bc_R[8] 0.640 0.150 0.357 0.640 0.941
mu_bc_R[9] 0.384 0.167 0.054 0.389 0.692
mu_bc_R[10] 0.892 0.102 0.690 0.892 1.092
mu_bc_R[11] 1.108 0.074 0.964 1.107 1.260
mu_bc_R[12] 0.889 0.160 0.568 0.893 1.202
mu_bc_R[13] 1.067 0.093 0.885 1.069 1.247
mu_bc_R[14] 0.854 0.121 0.613 0.855 1.079
mu_bc_R[15] 0.727 0.098 0.534 0.728 0.919
mu_bc_R[16] 1.270 0.116 1.037 1.270 1.496
tau_pH[1] 2.439 0.182 2.097 2.427 2.814
tau_pH[2] 1.585 0.148 1.314 1.582 1.881
tau_pH[3] 2.272 0.294 1.734 2.262 2.876
tau_pH[4] 1.053 0.146 0.791 1.043 1.363
tau_pH[5] 1.023 0.154 0.753 1.011 1.361
beta0_pH[1,1] 0.536 0.170 0.199 0.539 0.859
beta0_pH[2,1] 1.433 0.163 1.097 1.437 1.729
beta0_pH[3,1] 1.328 0.337 0.395 1.407 1.780
beta0_pH[4,1] 1.593 0.227 1.116 1.603 2.011
beta0_pH[5,1] -0.794 0.382 -1.735 -0.745 -0.179
beta0_pH[6,1] -0.591 0.474 -1.816 -0.461 0.020
beta0_pH[7,1] 0.752 0.010 0.732 0.752 0.773
beta0_pH[8,1] -0.749 0.492 -2.276 -0.663 -0.044
beta0_pH[9,1] -0.337 0.291 -1.017 -0.306 0.152
beta0_pH[10,1] 0.283 0.200 -0.121 0.287 0.658
beta0_pH[11,1] 0.069 0.326 -0.556 0.084 0.659
beta0_pH[12,1] 0.629 0.197 0.223 0.640 0.990
beta0_pH[13,1] -0.696 0.144 -0.987 -0.691 -0.429
beta0_pH[14,1] 0.025 0.171 -0.308 0.027 0.354
beta0_pH[15,1] 0.044 0.204 -0.444 0.068 0.378
beta0_pH[16,1] -0.899 0.209 -1.350 -0.888 -0.517
beta0_pH[1,2] 2.566 0.083 2.411 2.565 2.735
beta0_pH[2,2] 3.079 0.022 3.036 3.078 3.124
beta0_pH[3,2] 3.625 0.042 3.544 3.625 3.708
beta0_pH[4,2] 2.937 0.019 2.900 2.937 2.975
beta0_pH[5,2] 4.502 1.391 2.685 4.233 7.917
beta0_pH[6,2] 3.326 0.159 3.032 3.322 3.649
beta0_pH[7,2] 2.026 0.050 1.931 2.026 2.126
beta0_pH[8,2] 2.824 0.066 2.698 2.824 2.954
beta0_pH[9,2] 3.955 0.206 3.557 3.955 4.366
beta0_pH[10,2] 3.998 0.130 3.749 3.996 4.256
beta0_pH[11,2] 1.889 0.009 1.872 1.889 1.906
beta0_pH[12,2] 1.802 0.069 1.673 1.800 1.939
beta0_pH[13,2] 2.499 0.026 2.449 2.499 2.551
beta0_pH[14,2] 1.805 0.024 1.760 1.805 1.852
beta0_pH[15,2] 2.551 0.017 2.518 2.551 2.584
beta0_pH[16,2] 2.777 0.025 2.729 2.777 2.827
beta0_pH[1,3] 0.682 0.605 -1.013 0.844 1.367
beta0_pH[2,3] 2.126 0.017 2.093 2.127 2.160
beta0_pH[3,3] 2.689 0.024 2.642 2.688 2.736
beta0_pH[4,3] 2.980 0.034 2.915 2.980 3.046
beta0_pH[5,3] 3.874 2.856 0.993 3.047 11.164
beta0_pH[6,3] 1.254 0.333 0.668 1.235 1.968
beta0_pH[7,3] 0.857 0.049 0.760 0.855 0.954
beta0_pH[8,3] 0.285 0.120 0.032 0.291 0.501
beta0_pH[9,3] -2.866 1.106 -5.262 -2.740 -1.012
beta0_pH[10,3] -0.909 0.893 -2.990 -0.792 0.555
beta0_pH[11,3] 0.109 11.438 -22.747 0.272 22.626
beta0_pH[12,3] 0.155 11.580 -22.741 0.303 23.284
beta0_pH[13,3] 0.285 11.417 -22.140 0.483 22.926
beta0_pH[14,3] 0.079 11.489 -21.821 0.309 22.515
beta0_pH[15,3] 0.258 11.572 -21.861 0.369 23.250
beta0_pH[16,3] 0.216 11.349 -21.716 0.515 22.345
beta0_pH[11,4] 1.530 0.721 -0.295 1.618 2.720
beta0_pH[12,4] 1.283 1.088 -1.240 1.301 3.411
beta0_pH[13,4] 1.901 0.799 0.106 1.983 3.223
beta0_pH[14,4] 0.646 0.641 -0.909 0.721 1.707
beta0_pH[15,4] 0.838 0.563 -0.477 1.011 1.524
beta0_pH[16,4] 1.147 0.881 -1.139 1.212 2.628
beta0_pH[11,5] 1.288 0.105 1.102 1.277 1.521
beta0_pH[12,5] -2.318 0.423 -3.266 -2.279 -1.586
beta0_pH[13,5] 0.760 0.055 0.662 0.757 0.880
beta0_pH[14,5] 0.626 0.234 0.167 0.625 1.077
beta0_pH[15,5] 2.282 0.447 1.603 2.191 3.357
beta0_pH[16,5] -0.771 0.069 -0.903 -0.773 -0.635
beta1_pH[1,1] 2.561 0.448 1.828 2.516 3.564
beta1_pH[2,1] 1.587 0.332 1.008 1.565 2.365
beta1_pH[3,1] 2.002 0.669 1.130 1.799 3.748
beta1_pH[4,1] 1.903 0.481 1.183 1.828 3.249
beta1_pH[5,1] 2.305 0.492 1.517 2.236 3.494
beta1_pH[6,1] 2.314 0.942 1.129 2.027 4.636
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.542 1.035 1.165 2.219 4.786
beta1_pH[9,1] 2.134 0.751 1.219 1.946 4.436
beta1_pH[10,1] 2.122 0.370 1.491 2.088 2.928
beta1_pH[11,1] 3.365 0.540 2.472 3.296 4.614
beta1_pH[12,1] 2.905 0.292 2.347 2.906 3.518
beta1_pH[13,1] 3.517 0.230 3.079 3.506 3.997
beta1_pH[14,1] 2.833 0.309 2.228 2.828 3.457
beta1_pH[15,1] 2.415 0.289 1.907 2.392 3.072
beta1_pH[16,1] 4.079 0.288 3.544 4.065 4.691
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 3.435 1.639 1.230 2.949 7.237
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.865 2.322 2.202 4.101 11.261
beta1_pH[9,3] 8.439 2.631 4.186 8.104 14.038
beta1_pH[10,3] 7.246 2.742 3.195 6.767 13.375
beta1_pH[11,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,4] 1.981 0.765 0.566 1.929 3.732
beta1_pH[12,4] 4.117 1.117 2.013 4.078 6.672
beta1_pH[13,4] 2.487 0.829 1.051 2.408 4.294
beta1_pH[14,4] 2.763 0.700 1.582 2.706 4.387
beta1_pH[15,4] 0.911 0.603 0.220 0.735 2.314
beta1_pH[16,4] 3.907 0.929 2.299 3.825 6.303
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.290 0.842 0.358 0.976 3.343
beta2_pH[2,1] 1.674 0.958 0.308 1.493 3.432
beta2_pH[3,1] 1.227 0.907 0.139 0.997 2.942
beta2_pH[4,1] 1.144 0.793 0.209 0.868 2.908
beta2_pH[5,1] 2.255 1.895 0.293 1.756 7.795
beta2_pH[6,1] 1.041 1.741 0.103 0.373 6.606
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 0.859 1.543 0.108 0.324 6.277
beta2_pH[9,1] 1.098 1.574 0.133 0.543 6.062
beta2_pH[10,1] 1.448 1.686 0.303 0.864 7.167
beta2_pH[11,1] 0.425 0.237 0.181 0.364 1.053
beta2_pH[12,1] 1.530 0.798 0.610 1.270 3.647
beta2_pH[13,1] 0.942 0.352 0.534 0.860 1.849
beta2_pH[14,1] 1.099 0.632 0.460 0.915 3.115
beta2_pH[15,1] 1.027 0.626 0.361 0.853 2.914
beta2_pH[16,1] 1.028 0.306 0.626 0.975 1.786
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 2.209 4.080 0.120 0.498 15.548
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 2.549 2.194 0.571 1.591 8.503
beta2_pH[9,3] 0.190 0.190 0.081 0.158 0.453
beta2_pH[10,3] 0.303 0.347 0.102 0.226 0.957
beta2_pH[11,3] -8.537 6.063 -22.779 -7.410 -0.433
beta2_pH[12,3] -8.536 6.048 -22.504 -7.363 -0.482
beta2_pH[13,3] -8.542 6.112 -22.739 -7.282 -0.454
beta2_pH[14,3] -8.534 6.184 -22.643 -7.227 -0.427
beta2_pH[15,3] -8.530 6.072 -23.080 -7.333 -0.412
beta2_pH[16,3] -8.528 6.120 -22.883 -7.330 -0.473
beta2_pH[11,4] -4.494 3.616 -14.458 -3.423 -0.675
beta2_pH[12,4] -4.910 3.650 -14.651 -3.942 -0.887
beta2_pH[13,4] -4.215 3.633 -14.433 -3.023 -0.713
beta2_pH[14,4] -4.312 3.600 -14.666 -3.093 -0.806
beta2_pH[15,4] -4.688 3.596 -14.693 -3.633 -0.773
beta2_pH[16,4] -4.559 3.291 -13.878 -3.536 -1.167
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.784 1.059 32.825 34.779 37.086
beta3_pH[2,1] 32.783 1.372 30.514 32.731 35.836
beta3_pH[3,1] 33.326 2.035 29.292 33.146 38.599
beta3_pH[4,1] 32.836 2.035 29.674 32.623 37.758
beta3_pH[5,1] 27.558 1.086 25.935 27.455 30.302
beta3_pH[6,1] 33.394 3.243 26.719 33.179 40.649
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 34.642 4.064 26.955 34.537 42.171
beta3_pH[9,1] 32.411 3.304 27.793 31.705 41.249
beta3_pH[10,1] 32.535 1.301 30.253 32.439 35.323
beta3_pH[11,1] 31.379 1.103 29.322 31.420 33.552
beta3_pH[12,1] 30.803 0.484 29.813 30.821 31.715
beta3_pH[13,1] 31.382 0.557 30.359 31.342 32.558
beta3_pH[14,1] 32.558 0.750 31.211 32.527 34.131
beta3_pH[15,1] 30.906 0.729 29.391 30.941 32.273
beta3_pH[16,1] 30.122 0.366 29.389 30.116 30.833
beta3_pH[1,2] 34.393 5.532 25.484 34.411 43.556
beta3_pH[2,2] 34.482 5.462 25.520 34.450 43.610
beta3_pH[3,2] 30.400 7.180 18.577 30.453 42.321
beta3_pH[4,2] 30.390 7.217 18.591 30.501 42.331
beta3_pH[5,2] 30.366 7.138 18.485 30.389 42.334
beta3_pH[6,2] 30.459 7.193 18.726 30.402 42.443
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 39.653 2.713 33.083 40.141 43.677
beta3_pH[2,3] 34.327 5.482 25.403 34.356 43.570
beta3_pH[3,3] 30.328 7.173 18.639 30.251 42.457
beta3_pH[4,3] 30.281 7.108 18.575 30.335 42.369
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.120 0.823 40.932 41.871 43.861
beta3_pH[9,3] 35.928 1.640 34.527 35.289 40.577
beta3_pH[10,3] 37.777 2.497 34.629 37.113 43.338
beta3_pH[11,3] 38.414 2.586 34.205 38.375 42.728
beta3_pH[12,3] 38.482 2.558 34.268 38.470 42.730
beta3_pH[13,3] 38.576 2.604 34.257 38.524 42.814
beta3_pH[14,3] 38.449 2.620 34.208 38.409 42.751
beta3_pH[15,3] 38.464 2.628 34.241 38.451 42.793
beta3_pH[16,3] 38.481 2.586 34.223 38.375 42.761
beta3_pH[11,4] 40.568 1.206 38.493 40.396 43.244
beta3_pH[12,4] 42.247 0.670 40.456 42.392 43.099
beta3_pH[13,4] 40.973 0.977 39.170 40.952 42.985
beta3_pH[14,4] 40.243 0.733 38.854 40.219 41.759
beta3_pH[15,4] 41.306 1.434 37.784 41.441 43.766
beta3_pH[16,4] 41.739 0.470 40.877 41.769 42.711
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.263 0.021 0.221 0.262 0.305
beta4_pH[2,1] 0.211 0.005 0.201 0.211 0.221
beta4_pH[3,1] 0.172 0.007 0.159 0.172 0.186
beta4_pH[4,1] 0.294 0.011 0.272 0.294 0.316
beta4_pH[5,1] 0.424 0.301 0.083 0.385 0.994
beta4_pH[6,1] 0.789 0.392 0.316 0.740 1.489
beta4_pH[7,1] 0.848 0.342 0.336 0.794 1.650
beta4_pH[8,1] 0.156 0.069 0.042 0.149 0.310
beta4_pH[9,1] 0.215 0.086 0.062 0.211 0.395
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.447 0.012 0.424 0.446 0.470
beta4_pH[12,1] 0.643 0.023 0.598 0.642 0.689
beta4_pH[13,1] 0.629 0.013 0.603 0.629 0.656
beta4_pH[14,1] 0.669 0.031 0.612 0.668 0.730
beta4_pH[15,1] 0.796 0.021 0.755 0.796 0.839
beta4_pH[16,1] 1.214 0.041 1.136 1.213 1.299
beta4_pH[1,2] 0.296 0.096 0.124 0.289 0.502
beta4_pH[2,2] 0.514 0.023 0.470 0.514 0.559
beta4_pH[3,2] 0.446 0.012 0.423 0.446 0.469
beta4_pH[4,2] 0.654 0.038 0.579 0.653 0.731
beta4_pH[5,2] 0.736 0.989 0.115 0.539 2.537
beta4_pH[6,2] 0.910 0.848 0.219 0.736 2.647
beta4_pH[7,2] 0.897 0.538 0.349 0.811 1.874
beta4_pH[8,2] 0.850 0.202 0.475 0.842 1.273
beta4_pH[9,2] 1.150 1.039 0.286 0.883 3.758
beta4_pH[10,2] 0.457 0.052 0.359 0.455 0.563
beta4_pH[11,2] 1.270 0.050 1.176 1.269 1.371
beta4_pH[12,2] 0.611 0.111 0.410 0.609 0.838
beta4_pH[13,2] 0.953 0.045 0.866 0.952 1.045
beta4_pH[14,2] 0.935 0.064 0.813 0.933 1.062
beta4_pH[15,2] 1.009 0.025 0.960 1.009 1.061
beta4_pH[16,2] 1.391 0.077 1.246 1.389 1.545
beta4_pH[1,3] 0.241 0.078 0.102 0.236 0.404
beta4_pH[2,3] 0.050 0.017 0.020 0.049 0.085
beta4_pH[3,3] 0.052 0.010 0.034 0.052 0.072
beta4_pH[4,3] 0.167 0.031 0.110 0.167 0.229
beta4_pH[5,3] 0.747 0.895 0.117 0.547 2.566
beta4_pH[6,3] 0.880 0.783 0.208 0.725 2.360
beta4_pH[7,3] 0.896 0.431 0.357 0.816 1.938
beta4_pH[8,3] 0.390 0.182 0.107 0.367 0.817
beta4_pH[9,3] 0.918 0.858 0.166 0.710 2.797
beta4_pH[10,3] 0.055 0.038 0.003 0.048 0.138
beta4_pH[11,3] 0.629 0.280 0.254 0.584 1.299
beta4_pH[12,3] 0.730 0.250 0.379 0.687 1.322
beta4_pH[13,3] 0.801 0.380 0.313 0.722 1.761
beta4_pH[14,3] 0.861 0.220 0.505 0.835 1.382
beta4_pH[15,3] 0.797 0.204 0.470 0.774 1.259
beta4_pH[16,3] 0.828 0.373 0.324 0.764 1.685
beta4_pH[1,4] 0.635 0.402 0.201 0.553 1.541
beta4_pH[2,4] 0.830 0.907 0.113 0.592 2.957
beta4_pH[3,4] 0.846 0.890 0.117 0.603 3.111
beta4_pH[4,4] 0.727 0.457 0.202 0.623 1.858
beta4_pH[5,4] 0.774 0.993 0.119 0.562 2.721
beta4_pH[6,4] 0.908 0.727 0.229 0.734 2.859
beta4_pH[7,4] 0.892 0.428 0.367 0.818 1.871
beta4_pH[8,4] 0.739 0.650 0.154 0.595 2.335
beta4_pH[9,4] 0.742 0.740 0.120 0.549 2.565
beta4_pH[10,4] 1.424 2.849 0.006 0.382 10.232
beta4_pH[11,4] 0.427 0.019 0.391 0.427 0.464
beta4_pH[12,4] 0.534 0.045 0.445 0.533 0.622
beta4_pH[13,4] 0.836 0.033 0.774 0.836 0.900
beta4_pH[14,4] 0.877 0.061 0.765 0.875 1.003
beta4_pH[15,4] 0.617 0.018 0.582 0.617 0.653
beta4_pH[16,4] 0.665 0.044 0.580 0.665 0.755
beta4_pH[1,5] 0.633 0.377 0.196 0.558 1.520
beta4_pH[2,5] 0.850 1.033 0.111 0.584 3.233
beta4_pH[3,5] 0.853 0.893 0.115 0.594 3.242
beta4_pH[4,5] 0.737 0.507 0.208 0.620 1.912
beta4_pH[5,5] 0.745 0.895 0.108 0.552 2.613
beta4_pH[6,5] 0.923 0.805 0.203 0.751 2.575
beta4_pH[7,5] 0.886 0.427 0.358 0.804 1.844
beta4_pH[8,5] 0.736 0.690 0.149 0.587 2.089
beta4_pH[9,5] 0.793 0.984 0.119 0.561 2.854
beta4_pH[10,5] 1.587 3.345 0.008 0.397 12.685
beta4_pH[11,5] 0.311 0.081 0.161 0.308 0.477
beta4_pH[12,5] 0.865 0.286 0.466 0.806 1.592
beta4_pH[13,5] 2.550 0.352 1.870 2.543 3.274
beta4_pH[14,5] 0.846 0.195 0.514 0.830 1.280
beta4_pH[15,5] 0.922 0.075 0.795 0.914 1.097
beta4_pH[16,5] 0.964 0.398 0.387 0.898 1.950
beta0_pelagic[1] 2.267 0.036 2.197 2.266 2.340
beta0_pelagic[2] 1.538 0.007 1.524 1.538 1.552
beta0_pelagic[3] 0.207 0.231 -0.334 0.247 0.572
beta0_pelagic[4] 0.222 0.238 -0.422 0.241 0.612
beta0_pelagic[5] 1.441 0.154 1.150 1.439 1.754
beta0_pelagic[6] 1.355 0.061 1.238 1.353 1.481
beta0_pelagic[7] 1.470 0.037 1.401 1.468 1.546
beta0_pelagic[8] 1.672 0.043 1.594 1.671 1.761
beta0_pelagic[9] 2.774 0.063 2.659 2.771 2.901
beta0_pelagic[10] 2.603 0.031 2.544 2.602 2.666
beta0_pelagic[11] -0.137 0.180 -0.593 -0.116 0.148
beta0_pelagic[12] 1.799 0.035 1.733 1.798 1.872
beta0_pelagic[13] -2.266 1.194 -4.534 -2.141 -0.261
beta0_pelagic[14] -0.266 0.174 -0.662 -0.247 0.025
beta0_pelagic[15] -0.210 0.035 -0.304 -0.195 -0.180
beta0_pelagic[16] -0.940 0.634 -2.242 -0.864 -0.016
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.743 0.262 0.300 0.723 1.420
beta1_pelagic[4] 0.826 0.259 0.388 0.809 1.482
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.146 0.255 0.757 1.116 1.844
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 5.009 2.167 1.598 4.654 8.793
beta1_pelagic[14] 1.703 0.430 1.074 1.623 2.793
beta1_pelagic[15] 0.097 0.167 0.000 0.014 0.634
beta1_pelagic[16] 2.368 1.037 0.857 2.297 4.789
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 6.435 6.111 0.216 4.684 21.991
beta2_pelagic[4] 7.722 6.075 0.354 6.333 21.828
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 2.538 2.180 0.252 1.641 7.081
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.119 0.080 0.051 0.092 0.316
beta2_pelagic[14] 1.317 1.437 0.221 0.699 5.320
beta2_pelagic[15] 2.600 1.911 0.150 2.179 7.118
beta2_pelagic[16] 0.421 0.698 0.079 0.198 2.725
beta3_pelagic[1] 28.774 3.912 22.354 28.724 35.171
beta3_pelagic[2] 28.733 3.893 22.327 28.756 35.144
beta3_pelagic[3] 27.923 2.604 23.441 28.308 32.047
beta3_pelagic[4] 26.117 1.456 23.040 26.256 29.542
beta3_pelagic[5] 29.416 6.639 18.602 29.654 40.216
beta3_pelagic[6] 29.564 6.657 18.575 29.558 40.394
beta3_pelagic[7] 29.545 6.570 18.644 29.511 40.471
beta3_pelagic[8] 29.565 6.613 18.590 29.673 40.379
beta3_pelagic[9] 29.621 6.627 18.553 29.696 40.545
beta3_pelagic[10] 29.468 6.649 18.607 29.467 40.395
beta3_pelagic[11] 30.103 1.172 27.998 30.032 32.745
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 29.692 5.007 23.254 28.694 41.049
beta3_pelagic[14] 33.957 1.546 31.343 33.877 37.652
beta3_pelagic[15] 34.270 5.661 23.649 35.107 41.667
beta3_pelagic[16] 30.427 3.615 23.922 30.572 38.277
beta0_yellow[1] -0.483 0.080 -0.637 -0.482 -0.330
beta0_yellow[2] 0.539 0.087 0.339 0.549 0.680
beta0_yellow[3] -0.205 0.127 -0.460 -0.208 0.015
beta0_yellow[4] 0.910 0.234 0.253 0.973 1.191
beta0_yellow[5] 0.235 0.331 -0.476 0.252 0.826
beta0_yellow[6] 0.820 0.106 0.611 0.820 1.034
beta0_yellow[7] 1.010 0.102 0.812 1.010 1.216
beta0_yellow[8] 1.018 0.092 0.835 1.016 1.198
beta0_yellow[9] 0.442 0.099 0.250 0.442 0.634
beta0_yellow[10] 0.538 0.069 0.397 0.539 0.672
beta0_yellow[11] -0.225 0.229 -0.589 -0.186 0.002
beta0_yellow[12] -1.094 0.097 -1.281 -1.095 -0.908
beta0_yellow[13] -1.356 0.177 -1.732 -1.342 -1.065
beta0_yellow[14] -0.248 0.139 -0.565 -0.239 -0.014
beta0_yellow[15] -1.246 0.135 -1.574 -1.226 -1.035
beta0_yellow[16] -0.253 0.093 -0.428 -0.255 -0.072
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.019 0.238 0.592 1.007 1.495
beta1_yellow[3] 0.461 0.244 0.019 0.475 0.921
beta1_yellow[4] 1.330 0.581 0.659 1.185 3.160
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.257 0.350 0.001 0.074 1.172
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 0.999 1.052 0.019 0.725 4.756
beta1_yellow[14] 0.645 0.823 0.013 0.422 3.143
beta1_yellow[15] 0.166 0.219 0.011 0.080 0.740
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -7.672 5.867 -21.859 -6.370 -0.404
beta2_yellow[3] -7.719 5.861 -21.464 -6.438 -0.392
beta2_yellow[4] -5.596 5.831 -19.815 -3.644 -0.127
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -7.826 5.749 -21.424 -6.605 -0.423
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -6.149 5.155 -18.581 -5.059 -0.339
beta2_yellow[14] -5.908 4.975 -18.626 -4.432 -0.379
beta2_yellow[15] -7.686 5.720 -21.192 -6.575 -0.270
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.891 1.149 27.291 28.670 31.892
beta3_yellow[3] 32.328 3.356 23.138 32.503 38.832
beta3_yellow[4] 29.637 3.192 23.613 28.464 35.569
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 35.514 3.879 30.324 34.443 43.581
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 33.937 2.714 29.308 34.016 40.591
beta3_yellow[14] 33.633 3.389 29.175 33.201 42.526
beta3_yellow[15] 36.645 3.864 30.290 36.382 43.502
beta3_yellow[16] 35.000 0.000 35.000 35.000 35.000
beta0_black[1] 0.012 0.076 -0.139 0.012 0.159
beta0_black[2] 1.856 0.031 1.795 1.856 1.917
beta0_black[3] 1.239 0.037 1.165 1.239 1.311
beta0_black[4] 2.425 0.047 2.336 2.424 2.520
beta0_black[5] 5.232 2.496 1.728 4.722 11.416
beta0_black[6] 5.227 2.174 2.585 4.653 10.931
beta0_black[7] 4.104 2.138 1.633 3.517 10.055
beta0_black[8] 0.892 0.220 0.476 0.888 1.349
beta0_black[9] 2.818 0.266 2.315 2.818 3.362
beta0_black[10] 1.288 0.050 1.194 1.287 1.389
beta0_black[11] 2.978 0.052 2.875 2.977 3.080
beta0_black[12] 4.908 0.196 4.519 4.904 5.297
beta0_black[13] 0.446 0.340 0.016 0.366 1.224
beta0_black[14] 2.852 0.103 2.653 2.852 3.058
beta0_black[15] 0.722 0.052 0.621 0.722 0.824
beta0_black[16] 3.772 0.073 3.629 3.771 3.912
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.525 1.763 -6.695 -2.136 -0.287
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.265 3.482 26.132 38.109 41.013
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.327 0.092 -0.511 -0.324 -0.147
beta4_black[2] 0.258 0.052 0.158 0.257 0.363
beta4_black[3] -0.718 0.064 -0.849 -0.718 -0.594
beta4_black[4] 0.632 0.130 0.390 0.630 0.900
beta4_black[5] 0.846 1.862 -1.085 0.531 4.613
beta4_black[6] 0.927 1.819 -0.909 0.576 4.684
beta4_black[7] 0.857 1.800 -1.040 0.523 4.854
beta4_black[8] 0.081 0.287 -0.526 0.091 0.611
beta4_black[9] 1.063 0.837 -0.103 0.884 3.201
beta4_black[10] 0.074 0.066 -0.057 0.075 0.206
beta4_black[11] -0.611 0.065 -0.734 -0.613 -0.481
beta4_black[12] 1.232 0.290 0.679 1.220 1.813
beta4_black[13] -1.003 0.077 -1.161 -1.001 -0.856
beta4_black[14] 0.216 0.180 -0.125 0.210 0.582
beta4_black[15] 0.003 0.046 -0.088 0.003 0.094
beta4_black[16] -0.416 0.115 -0.641 -0.418 -0.181
beta0_dsr[11] -0.338 0.009 -0.356 -0.338 -0.320
beta0_dsr[12] 1.420 0.044 1.336 1.420 1.508
beta0_dsr[13] 0.143 0.022 0.100 0.143 0.187
beta0_dsr[14] -0.129 0.017 -0.163 -0.128 -0.095
beta0_dsr[15] 0.196 0.014 0.169 0.196 0.223
beta0_dsr[16] 0.208 0.015 0.178 0.208 0.237
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.548 0.035 0.480 0.548 0.616
beta4_dsr[12] 0.796 0.233 0.355 0.789 1.257
beta4_dsr[13] -0.078 0.040 -0.154 -0.079 0.001
beta4_dsr[14] -0.213 0.072 -0.361 -0.211 -0.076
beta4_dsr[15] 0.225 0.021 0.183 0.225 0.268
beta4_dsr[16] -0.099 0.046 -0.188 -0.098 -0.010
beta0_slope[11] -2.197 0.021 -2.239 -2.197 -2.157
beta0_slope[12] -5.156 0.172 -5.474 -5.163 -4.810
beta0_slope[13] -1.474 0.033 -1.540 -1.473 -1.409
beta0_slope[14] -2.683 0.047 -2.777 -2.683 -2.590
beta0_slope[15] -1.591 0.022 -1.635 -1.590 -1.547
beta0_slope[16] -3.016 0.046 -3.109 -3.015 -2.930
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.345 0.076 -0.494 -0.343 -0.200
beta4_slope[12] -4.246 0.815 -5.950 -4.202 -2.836
beta4_slope[13] 0.308 0.053 0.204 0.309 0.413
beta4_slope[14] -0.192 0.152 -0.499 -0.188 0.097
beta4_slope[15] -0.363 0.033 -0.431 -0.364 -0.299
beta4_slope[16] 0.221 0.104 0.014 0.223 0.418
beta5_slope[11] 3.071 0.106 2.874 3.068 3.285
beta5_slope[12] 1.578 0.403 0.818 1.554 2.397
beta5_slope[13] 0.847 0.077 0.704 0.846 0.999
beta5_slope[14] 3.785 0.357 3.117 3.771 4.541
beta5_slope[15] 1.416 0.073 1.275 1.415 1.563
beta5_slope[16] 2.684 0.233 2.255 2.676 3.163
tau_comp[1] 3.781 0.298 3.236 3.763 4.384
tau_comp[2] 2.858 0.257 2.405 2.845 3.391
tau_comp[3] 0.869 0.075 0.730 0.868 1.023
tau_comp[4] 1.358 0.214 0.984 1.344 1.816
tau_comp[5] 1.205 0.178 0.887 1.194 1.578
sigma_H[1] 0.207 0.047 0.125 0.203 0.306
sigma_H[2] 0.175 0.030 0.125 0.172 0.241
sigma_H[3] 0.190 0.040 0.120 0.188 0.273
sigma_H[4] 0.460 0.080 0.329 0.451 0.645
sigma_H[5] 0.962 0.204 0.588 0.950 1.381
sigma_H[6] 0.538 0.195 0.160 0.535 0.935
sigma_H[7] 0.310 0.062 0.216 0.302 0.453
sigma_H[8] 0.419 0.089 0.274 0.408 0.619
sigma_H[9] 0.514 0.119 0.332 0.501 0.784
sigma_H[10] 0.211 0.039 0.143 0.207 0.297
sigma_H[11] 0.265 0.044 0.192 0.260 0.364
sigma_H[12] 0.382 0.145 0.197 0.339 0.713
sigma_H[13] 0.220 0.036 0.160 0.216 0.299
sigma_H[14] 0.499 0.089 0.348 0.491 0.691
sigma_H[15] 0.229 0.037 0.166 0.226 0.310
sigma_H[16] 0.213 0.041 0.146 0.210 0.303
lambda_H[1] 2.338 3.631 0.080 1.033 12.263
lambda_H[2] 5.088 4.900 0.427 3.645 17.464
lambda_H[3] 5.268 7.416 0.258 2.784 24.476
lambda_H[4] 0.007 0.005 0.001 0.005 0.018
lambda_H[5] 2.022 4.190 0.027 0.625 13.552
lambda_H[6] 6.745 14.282 0.010 0.499 47.261
lambda_H[7] 0.012 0.009 0.002 0.010 0.034
lambda_H[8] 7.956 9.811 0.102 4.403 34.589
lambda_H[9] 0.015 0.010 0.003 0.013 0.040
lambda_H[10] 0.242 0.493 0.027 0.162 0.863
lambda_H[11] 0.219 0.330 0.010 0.101 1.052
lambda_H[12] 3.629 4.752 0.194 2.048 16.765
lambda_H[13] 3.241 2.895 0.271 2.438 10.880
lambda_H[14] 3.342 4.201 0.223 2.067 14.584
lambda_H[15] 0.019 0.018 0.003 0.014 0.062
lambda_H[16] 0.578 0.869 0.034 0.299 2.698
mu_lambda_H[1] 4.086 1.861 0.996 3.927 8.025
mu_lambda_H[2] 3.690 1.961 0.500 3.511 7.890
mu_lambda_H[3] 3.329 1.833 0.711 3.044 7.688
sigma_lambda_H[1] 8.167 4.401 1.590 7.464 17.968
sigma_lambda_H[2] 8.136 4.737 0.841 7.529 18.493
sigma_lambda_H[3] 6.122 4.064 0.917 5.140 16.446
beta_H[1,1] 6.538 1.307 3.504 6.789 8.376
beta_H[2,1] 9.759 0.630 8.368 9.784 10.925
beta_H[3,1] 7.822 0.801 6.014 7.907 9.134
beta_H[4,1] 10.746 7.620 -4.962 10.976 25.438
beta_H[5,1] -0.447 2.642 -5.726 -0.307 4.093
beta_H[6,1] 3.271 4.193 -7.166 4.804 8.069
beta_H[7,1] 1.406 5.633 -10.745 1.716 11.466
beta_H[8,1] 1.280 3.440 -2.543 1.209 3.557
beta_H[9,1] 13.623 5.644 2.363 13.549 24.989
beta_H[10,1] 6.881 1.782 3.116 6.929 10.212
beta_H[11,1] 4.780 3.573 -2.929 5.374 9.978
beta_H[12,1] 2.624 1.107 0.697 2.541 5.083
beta_H[13,1] 9.070 0.874 7.042 9.137 10.559
beta_H[14,1] 2.121 1.033 0.010 2.150 4.062
beta_H[15,1] -7.037 3.361 -13.240 -7.077 -0.005
beta_H[16,1] 3.916 2.965 -0.862 3.498 10.841
beta_H[1,2] 7.781 0.272 7.231 7.786 8.284
beta_H[2,2] 9.884 0.163 9.560 9.883 10.203
beta_H[3,2] 8.798 0.201 8.410 8.799 9.201
beta_H[4,2] 2.992 1.458 0.163 2.944 5.891
beta_H[5,2] 1.685 0.997 -0.278 1.686 3.601
beta_H[6,2] 6.110 1.049 3.604 6.263 7.743
beta_H[7,2] 2.358 1.034 0.487 2.287 4.566
beta_H[8,2] 3.005 1.028 1.316 3.141 4.237
beta_H[9,2] 3.221 1.078 1.167 3.201 5.412
beta_H[10,2] 8.037 0.362 7.291 8.045 8.724
beta_H[11,2] 9.844 0.649 8.856 9.741 11.252
beta_H[12,2] 3.932 0.366 3.237 3.921 4.698
beta_H[13,2] 9.146 0.236 8.707 9.139 9.626
beta_H[14,2] 4.028 0.345 3.344 4.033 4.688
beta_H[15,2] 11.538 0.616 10.294 11.556 12.691
beta_H[16,2] 4.519 0.815 2.989 4.490 6.132
beta_H[1,3] 8.431 0.260 7.961 8.426 8.939
beta_H[2,3] 9.858 0.147 9.577 9.854 10.153
beta_H[3,3] 9.460 0.160 9.151 9.457 9.788
beta_H[4,3] -2.129 0.866 -3.888 -2.120 -0.521
beta_H[5,3] 3.870 0.628 2.592 3.882 5.043
beta_H[6,3] 8.289 1.135 6.690 7.957 10.657
beta_H[7,3] -2.540 0.614 -3.760 -2.546 -1.321
beta_H[8,3] 5.274 0.498 4.665 5.204 6.290
beta_H[9,3] -2.752 0.724 -4.232 -2.727 -1.363
beta_H[10,3] 8.542 0.274 8.019 8.534 9.092
beta_H[11,3] 8.495 0.293 7.873 8.521 9.016
beta_H[12,3] 5.269 0.327 4.503 5.305 5.817
beta_H[13,3] 8.876 0.175 8.516 8.877 9.216
beta_H[14,3] 5.781 0.270 5.160 5.804 6.251
beta_H[15,3] 10.271 0.299 9.685 10.264 10.855
beta_H[16,3] 6.045 0.585 4.822 6.077 7.063
beta_H[1,4] 8.093 0.187 7.697 8.101 8.435
beta_H[2,4] 9.847 0.130 9.588 9.852 10.095
beta_H[3,4] 9.951 0.151 9.616 9.964 10.220
beta_H[4,4] 11.806 0.468 10.866 11.820 12.704
beta_H[5,4] 5.591 0.736 4.391 5.515 7.237
beta_H[6,4] 7.246 0.908 5.297 7.499 8.498
beta_H[7,4] 8.536 0.345 7.859 8.545 9.186
beta_H[8,4] 6.723 0.242 6.224 6.736 7.126
beta_H[9,4] 7.440 0.454 6.558 7.438 8.338
beta_H[10,4] 7.536 0.226 7.101 7.528 8.014
beta_H[11,4] 9.502 0.191 9.139 9.499 9.878
beta_H[12,4] 7.245 0.195 6.870 7.238 7.652
beta_H[13,4] 9.062 0.141 8.772 9.068 9.329
beta_H[14,4] 7.844 0.218 7.439 7.840 8.277
beta_H[15,4] 9.461 0.221 9.029 9.461 9.890
beta_H[16,4] 9.399 0.237 8.938 9.392 9.886
beta_H[1,5] 8.894 0.152 8.585 8.894 9.194
beta_H[2,5] 10.674 0.100 10.487 10.672 10.877
beta_H[3,5] 10.762 0.162 10.465 10.758 11.098
beta_H[4,5] 8.361 0.507 7.400 8.352 9.377
beta_H[5,5] 5.181 0.617 3.719 5.258 6.166
beta_H[6,5] 8.938 0.627 8.002 8.798 10.314
beta_H[7,5] 6.575 0.343 5.914 6.570 7.276
beta_H[8,5] 8.255 0.206 7.909 8.244 8.705
beta_H[9,5] 8.202 0.456 7.280 8.203 9.090
beta_H[10,5] 10.034 0.223 9.593 10.036 10.467
beta_H[11,5] 11.596 0.219 11.161 11.594 12.023
beta_H[12,5] 8.614 0.191 8.242 8.610 9.003
beta_H[13,5] 10.127 0.132 9.875 10.126 10.398
beta_H[14,5] 9.245 0.216 8.850 9.233 9.702
beta_H[15,5] 11.301 0.235 10.843 11.299 11.762
beta_H[16,5] 9.875 0.184 9.496 9.882 10.222
beta_H[1,6] 10.123 0.239 9.756 10.086 10.703
beta_H[2,6] 11.474 0.115 11.246 11.473 11.715
beta_H[3,6] 10.623 0.156 10.287 10.630 10.912
beta_H[4,6] 12.898 0.899 11.062 12.917 14.606
beta_H[5,6] 5.660 0.623 4.556 5.631 7.027
beta_H[6,6] 8.746 0.548 7.511 8.781 9.753
beta_H[7,6] 10.547 0.589 9.374 10.540 11.680
beta_H[8,6] 9.536 0.255 9.040 9.551 9.946
beta_H[9,6] 8.268 0.763 6.801 8.255 9.760
beta_H[10,6] 9.154 0.340 8.413 9.178 9.756
beta_H[11,6] 10.673 0.359 9.887 10.697 11.285
beta_H[12,6] 9.369 0.236 8.927 9.365 9.870
beta_H[13,6] 11.071 0.159 10.771 11.067 11.389
beta_H[14,6] 9.921 0.286 9.333 9.923 10.465
beta_H[15,6] 10.522 0.423 9.689 10.523 11.351
beta_H[16,6] 10.583 0.253 10.037 10.597 11.061
beta_H[1,7] 10.143 1.239 6.939 10.441 11.804
beta_H[2,7] 12.113 0.535 10.966 12.146 13.074
beta_H[3,7] 10.295 0.668 8.789 10.364 11.432
beta_H[4,7] 1.536 4.552 -6.886 1.481 10.792
beta_H[5,7] 5.845 2.019 1.971 5.893 9.727
beta_H[6,7] 7.377 3.042 -0.587 8.470 11.119
beta_H[7,7] 7.646 2.948 1.926 7.617 13.610
beta_H[8,7] 10.872 0.900 9.189 10.880 12.529
beta_H[9,7] 5.090 3.949 -2.715 5.130 12.878
beta_H[10,7] 10.019 1.612 7.101 9.885 13.596
beta_H[11,7] 11.401 1.829 8.273 11.187 15.516
beta_H[12,7] 9.763 0.972 7.604 9.824 11.462
beta_H[13,7] 11.668 0.731 10.004 11.740 12.835
beta_H[14,7] 10.634 0.917 8.697 10.664 12.388
beta_H[15,7] 14.067 2.201 9.819 14.079 18.519
beta_H[16,7] 12.815 1.334 10.561 12.695 15.737
beta0_H[1] 8.275 18.210 -27.671 8.637 45.170
beta0_H[2] 10.428 8.789 -6.962 10.477 27.248
beta0_H[3] 10.191 10.265 -10.200 9.962 32.422
beta0_H[4] 7.631 183.675 -362.516 5.898 376.413
beta0_H[5] 3.955 29.846 -51.555 4.283 60.087
beta0_H[6] 7.025 50.796 -100.936 7.797 111.248
beta0_H[7] 4.255 142.648 -278.289 4.691 293.210
beta0_H[8] 6.702 30.360 -18.355 6.447 28.763
beta0_H[9] 8.069 120.931 -238.223 8.121 247.088
beta0_H[10] 8.881 36.673 -69.439 9.002 85.040
beta0_H[11] 10.256 52.138 -105.076 10.200 125.216
beta0_H[12] 6.292 12.040 -18.531 6.621 30.525
beta0_H[13] 9.563 11.378 -11.861 9.907 30.375
beta0_H[14] 7.400 11.525 -15.085 7.206 32.889
beta0_H[15] 5.654 116.114 -225.201 5.473 247.925
beta0_H[16] 9.451 31.318 -53.831 8.480 73.300